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A physical model is proposed for the estimation of the screech amplitude in 
underexpanded supersonic jets. The model is based on the hypothesis that the 
interaction of a plane acoustic wave with stationary shock waves provides 
amplification of the transmitted acoustic wave upon traversing the shock. 
Powell’s discrete source model for screech incorporating a stationary array of 
acoustic monopoles is extended to accommodate variable source strength. The 
proposed model reveals that the acoustic sources are of increasing strength with 
downstream distance. It is shown that the screech amplitude increases with the 
fully expanded jet Mach number. Comparisons of predicted screech amplitude 
with available test data show satisfactory agreement. The effect of variable source . 
strength on directivity of the fundamental (first harmonic, lowest frequency 
mode) and the second harmonic (overtone) is found to be unimportant with regard 
to the principal lobe (main or major lobe) of considerable relative strength, and is 
appreciable only in the secondary or minor lobes (of relatively weaker strength). 

INTRODUCTION 

Supersonic jet screech constitutes an important technological challenge in the theory of jet 
noise. Screech refers to high-amplitude discrete tones superposed on the otherwise broadband 
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acoustic spectrum observed in imperfectly expanded supersonic jets (Fig. 1). The shock noise 
produced by imperfectly expanded (choked) jets is very complex, and consists of two separate 
components, namely broadband shock noise (arising from shock-turbulence interaction) and 
screech tones. This shock noise is in addition to the broadband turbulent mixing noise in the jet 
(Lighthill 1952, 1954; Ffowcs Williams 1963). In the nearfield, the very high intensity of these 
screech tones (as high as 170 to 180 dB) with a highly directional character with significant 
upstream propagating component can induce fatigue and cause serious structural damage and 
failure of aerodynamic vehicles in flight. A fundamental understanding of jet screech amplitude 
is thus requisite to the design of a new generation of rocket engines (Kandula 2006a, b), turbo-jet 
aircraft engines, and to a number of other flow circumstances such as jet impingement and cavity 
resonance involving feedback cycles. 

Owing to pressure mismatch (nozzle exit pressure exceeds the ambient pressure), 
underexpanded supersonic jet structures manifest themselves in shock-cell structures involving 
normal/oblique shocks and regions of expansion (Fig. 2). Under certain conditions, a Mach disc 
is also formed. These shock-expansion units interact with instability waves, vortices, turbulence, 
and other stream disturbances in the viscous shear layer that surrounds the inviscid region. 

Powell (1953a-c) first identified the mechanism of supersonic jet screech in terms of a 
resonant feedback loop or cycle that is self-sustaining and consists of the following key 
processes (Fig. 2): downstream passage of flow disturbances in the jet shear layer initiated near 
the nozzle lip region (where the relatively thin shear layer is susceptible to instability); 
interaction of these stream disturbances with the nearly periodic and stationary oblique shock- 
cell structure of the jet resulting in their amplification, with the ultimate generation and emission 
of intense acoustic waves or energy; upstream propagation of the acoustic waves immediately 


2 



outside the jet shear layer (acoustic feedback); excitation and maintenance of stream disturbances 
in the vicinity of the nozzle lip by the upstream propagating acoustic waves (so-called 
receptivity), which then propagate downstream, thus closing the feedback cycle without any 
external interference. Experimental data (for example, flow visualization measurements by Kaji 
& Nishijima 1996) suggest that the shock-noise is most intense somewhere between the third and 
fourth shock-cell structures. For a detailed discussion on screech, Seiner (1984) may be 
consulted. 

Powell’s theory of stationary array of discrete monopole sources distributed through the 
shock-cells (Powell 1953a) provided directivity patterns of screech that agree with 
measurements, as reported by Norum (1983) for unheated jets. It also furnished a framework for 
the estimation of screech frequency based on shock-cell spacing. Fisher & Morfey (1976) and 
Tam et al. (1986) provided expressions for the prediction of screech frequency over a wide range 
of fully expanded jet Mach number M y and jet total temperature T 0 . Tam’s expression for the 

screech is predicated on the so-called weakest-link theory, according to which the feedback 
acoustic waves represent the weakest link and control the screech frequency. Recently, Massey 
& Ahuja(1997) expanded on Tam’s formulation and provided expressions for screech frequency 
for hot jets with different screech modes. The screech mode changes from axisymmetric to 
helical and flapping modes as the Mach number increases. 

On the basis of Powell’s model of stationary point sources, Harper-Bourne and Fisher (1973) 
experimentally and theoretically treated broadband shock noise (radiated to the farfield) by 
considering shock-turbulence interaction. Harper-Bourne and Fisher considered spatial 
coherence of the disturbance field to account for the broadband sound radiation to the farfield. 
Generally speaking, the intensity of broadband shock-associated noise is more prevalent in the 
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forward arc, and is a function only of the pressure ratio. A stochastic model for broadband shock 
noise was reported by Tam (1987) who considered the nonlinear interaction between jet 
turbulence and stationary shock-cell system. Shock-turbulence interaction (shear layer and an 
oblique compression wave) was numerically investigated by Lui & Lele (2002) including the 
dynamic motion of the shock during interaction with shear layer vortices. 

With regard to the screech process, our current understanding of the highly nonlinear flow 
dynamics and of the rate of amplification of stream disturbances is rather incomplete. Whereas 
screech frequency is controlled by the feedback acoustic waves, the screech intensity is governed 
primarily by the characteristics of the downstream propagating flow disturbances and their 
interaction with the shock-cell structures. Several detailed measurements (Kaji & Nishijima 
1996; Norum 1982; Panda et al. 1997; Panda 1998, 1999; Massey et al. 1993; Krothapalli et al. 
1986; Petitjean 2005) and numerical viscous simulations (Kaji & Nishijima 1996; Shen & Tam 
2002; Manning & Lele 1998, 2002; Loh et al. 2003; Li & Gao 2004, 2005) have been reported 
and served to throw light on the subject. The experimental description of A-D mode switch 
(A l ,A 2 ,B,C,D; Norum 1982) according to the jet operating parameters represents a notable step. 
Despite the theoretical advances in screech frequency predictions, directivity character, 
modulations of instability waves by shocks (Glaznev 1973; Sedel’nikhov 1967) and unsteady 
shock motions and their role in noise generation, screech amplitude predictions remain elusive 
(1998, 1999). Owing to strong nonlinearity of the feedback loop, even an empirical formula for 
screech intensity is unavailable (Shen & Tam 2002). 

The preponderating questions concerning the interaction of shock waves with stream 
disturbances and the excitation of disturbances by upstream-traveling acoustic waves largely 
remain unanswered. The purpose of the present report is to investigate the amplitude of screech 
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based on a physical model for the interaction of stream disturbances (primarily regarded as 
finite-amplitude acoustic waves) with shock waves. This paper is primarily based on the author’s 
recent work (Kandula 2006c). 

L PROPOSED PHYSICAL MODEL 

A. Source Distribution 

The present model builds on Powell’s conception of stationary discrete acoustic monopoles 
(Fig. 3). However, in the present context account is taken of the variable strength of the 
monopoles, as the flow disturbances are amplified on their passage through the shock-cell 
system. As with Harper-Bourne and Fisher (1973) the acoustic sources are regarded to be 
distributed along the lip line at locations in the jet shear layer where shock reflections take place. 
It is assumed that the sources S U S 2 ,S 3 are equally spaced. It is also assumed that relative to the 
source S 2 , the source 5, leads in phase r and S 3 lags in an equal amount of* phase r . 

The downstream-propagating stream disturbances near the supersonic side of the shear layer 
and the upstream-traveling acoustic disturbances on the subsonic outer edge of the shear layer 
form two key components of the feedback process. Attention is focused here on the amplification 
of downstream-traveling stream disturbances. 

B. Shock-Acoustic Wave Interaction 

1. Transmission of an Acoustic Wave Through a Normal Shock 

The proposed model for the screech amplitude is based on the doctrine that the downstream- 
propagating stream disturbances in the screech cycle are regarded as acoustic-like (nearly 
isentropic pressure fluctuations) and are of finite amplitude. The amplification of these 
disturbances during their passage is a consequence of the interaction between an acoustic wave 
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(of finite amplitude) and a stationary shock wave. If the screech cycle is regarded as analogous to 
a thermodynamic cycle (with energy transfers occasioned during shock interaction and scattering 
at the nozzle lip), the acoustic waves form the working medium for the screech cycle. Thus, it is 
likely that the upstream-propagating acoustic waves and downstream-traveling acoustic waves 
constitute the screech cycle. The scattering of upstream-propagating acoustic waves by the 
nozzle lip ensures that this is so, even though some of the incoming acoustic energy is 
transmitted to the hydrodynamic disturbances or instability waves in the shear layer in the 
receptivity process. 

Unsteady phase-averaged measurements of nearfield pressure fluctuations in a round jet at 
M j =1.19 by Panda (1996) indeed reveal that the downstream-propagating disturbances contain 

both hydrodynamic pulsations and acoustic waves, whereas upstream-propagating components 
are only acoustic (Raman 1999). 

When an acoustic wave (traversing with sound speed relative to the moving fluid) is incident 
upon the shock from ahead of it (Fig. 4), no reflected acoustic wave can arise because the flow 
ahead is supersonic, and the acoustic wave (longitudinal isentropic sound wave) is transmitted 
(refracted) downstream of the shock. At the same time, perturbations in entropy and vorticity 
(both are nonisentropic and nonacoustic; these isobaric, transverse waves are termed entropy- 
vorticity waves after Carrier) are generated, which are transported along with the flow behind the 
shock (Mckenzie & Westphal 1968; Kontorovich 1960). On the other hand, when an acoustic 
wave is incident upon the shock from behind it, no disturbances can arise in the medium ahead of 
the shock. The acoustic waves are merely reflected, and the entropy and vorticity waves are 
generated behind the shock. Thus, the refraction of the incident acoustic wave is assumed as long 
as the approach flow is supersonic. 


6 



The main relationships developed here are based on transmission through normal shocks, 
although the shocks in the jet are not wholly normal but oblique. The Mach number normal to 
the shock M\ n is related to the approach Mach number M\ by the relation 

M ]„ = M\ sin /?] 

where (3\ represents the shock wave angle (angle between the approach velocity and the plane of 
the shock). However, it is believed that the relationships based on the normal-shock assumption 
provide an upper limit to the screech tone intensity. 

The passage and linear interaction of plane disturbances (sound waves) with a shock wave 
upon normal incidence for an ideal gas was first treated by Blokhitsev (1945) and Burgers 
(1946). When the perturbing waves are normally incident upon the shock, the transmission 
coefficient is expressed by Blokhitsev (1945) and Burgers (1946) as 

(la) 

where 5p x and Sp 2 denote the strength of the incident and refracted acoustic waves, 
respectively, M ln the normal Mach number of the flow upstream of the shock, M 2n the normal 
Mach number of the flow downstream of the shock, and y the specific heat ratio. The 
downstream Mach number is obtained from the Rankine-Hugonoit normal shock relation 
(Shapiro 1953) 
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Fig. 5 presents the variation of the pressure amplification Sp 2 / Sp\ as a function of incoming 
Mach number Mj. The pressure amplification is seen to rapidly grow as the upstream Mach 

number is increased. For large values ofMj , the pressure ratio varies asMj . For oblique shocks 
(/? < 90 deg), the pressure amplification is less than that for the normal shock. 

The result of the linear shock-acoustic wave interaction theory has been shown to be in close 
agreement with the numerical nonlinear Euler simulations by Zang et a. (1984) for incident 
angles within about 30 deg of the critical angle (zero degree incident angle refers to normal 
incidence) for the acoustic as well as vorticity wave. Although linear theory is employed here 
shock-acoustic wave interaction, the successive amplification is considered nonlinear. 

2. Successive Amplification (Shock-Cell Interaction) 

Based on the general behavior of the Mach number decay for typical underexpanded jet 
conditions, the following distribution of the approach (or upstream) Mach number for each tip of 
the shock-expansion unit (shock-cell reflection occurring in the shear layer) may be considered 
as a rough approximation: 

My = M r ~{M r - 1 / sin /?] )(i - l)/(m - 1) i = \,2,...m (2a) 

where M r is the upstream Mach number for the first shock-cell, My the incident Mach number 
in the i th shock-cell, and m the total number of active shock-cells under consideration. Thus, 
Eq. (2a) ensures that My ~ M r for the first shock-cell (/ = 1 ), and My = 1 /sin fl\ for the last 

cell ( / = m). Based on the available data, the quantity M r is approximated as 

M r « 1 AM j (2b) 


8 



The quantity My could be more accurately obtained by a direct calculation, using numerical 
methods. In the present work, the choice of Eq. (2a) for approximating My is based on the 
consideration that it facilitates a simpler formulation of the model. 

The fully expanded Mach number M is related to the pressure ratio p t / p a and the 

pressure ratio parameter /3 by 



fp \ir-tyr 

\Pa j 


-1 


(2c) 


and 



(2d) 


with p t and p a denoting jet total pressure and the ambient pressure, respectively. The Prandtl- 
Glauert parameter (3 is thus related to the strength of the normal shock wave. The number of 
shock reflections for a convergent-divergent nozzle strongly depends on the nozzle design Mach 
number M d and the operating pressure ratio. 

The relative source strength in each shock-cell can be determined from Eqs. (1) and (2). Eqs. 
(la) and (2a) assert that there is successive amplification of the sound waves during their passage 
through the shock-cell system, with the highest amplification rate occasioned in the first shock- 
cell. The overall (or total) amplification becomes 
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where A p in refers to the disturbance level upstream of the first shock-cell in the vicinity of the 
nozzle exit. 

II. SCREECH CHARACTER 

A. Screech Amplitude 

The screech amplitude (so-called excess tone level) relative to the broadband spectrum SPL s 
is expressed as 

SPL s = 20 log! o {Ap m / A p in ) (4) 

where A p m / Ap in is determined from Eq. (3). Thus, Eqs. (la) to (4) enable us to compute the 

excess (or relative) screech amplitude as a function of jet Mach number. The excess tone noise 
represents the difference between the screech amplitude and the broadband level (both shock- 
induced and turbulence-induced) at the screech frequency in a given spectrum (Massey et al. 
1993). 

The above theory does not display a separate effect of jet temperature on screech tone 
intensity. Data of Massey et al (1993) suggest that the tone intensity reduces with an increase in 
temperature. According to the measurements of Massey and Ahuja (1997) mode switching 
occurs at elevated jet temperature, thereby affecting screech intensity. 

B. Screech Directivity 

In the present formulation, an extension of Powell’s treatment (Powell 1953a, 1953b, 1953c) 
is considered for the directivity of screech based on stationary point sources along a line 
(Rayleigh 1888). The farfield radiation is determined synthetically by combining these sources. 
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With each source of strength of 4 n and 4 nH for the fundamental and the first harmonic 
respectively, the velocity potential distribution for a point source is of the form 

< fi = }- e Hct-r) + H_ e ik(ct-r) (5) 

r r 

where c denotes the sound speed, r the radial distance, t the time elapsed, and k the wave 
number, which is related to the wavelength A of the feedback acoustic wave by 

A = 2n / k (6) 

Following Powell, we assume that the phases rj = r 2 = r. For small ratios of sir (where 
s denotes the shock-cell spacing), the potential for the fundamental of a three-cell shock system 
(for example) can be expressed as 

(f)j- - ^ + ° 2 — — ^ D y cos k(ct - r-S t ) (7) 

The quantities «2 an ^ a 3 respectively refer to the relative strengths of the second and the third 

shock-cell with the strength of the first shock-cell taken as unity. The phasing between 
consecutive monopoles is determined by shock-cell spacing and disturbance convection Mach 
number M c (Norum 1982). 

The directivity D f can be expressed as 



- \ +(l + 03 )cos&(cr-.scos#)] 

(l + a 2 +a 3 ) 

(8) 

with 

T = s/u c - S /(cM c ) 

(9) 
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where 9 denotes the angle from the downstream flow direction and u c and M c refer to 
convective velocity and convective Mach number of the disturbance, respectively. Eq. (8) can 
be recast as 


; (l + «2 +a 3) 


2 ^ 

a 2 +(l + fl3)cos (l — M c cos 6) 

XM r 


( 10 ) 


Thus, the directionality of screech strongly depends on the shock cell spacing and the 
disturbance Mach number (Powell 1953a). The quantity (\-M c cos 9) is a Doppler factor 
incorporating the variation in retarded time and source phasing (Harper-Bourne & Fisher (1973). 
An expression for the directivity of the harmonic D h is provided by Eq. (10) with 2 n replaced 

by 4 n . In the special case ofa2 = a 3 =1, Eq. (10) reduces to that derived by Powell (1953a) for 

uniform source strength. Eq. (10) can be extended to four or more sources. 

The shock-cell spacing s is obtained from experimental data and is related to the pressure 
ratio parameter /? defined by 

s/dj=f<J3) 


For example, data of Harper-Bourne and Fisher (1973) suggest that 

s/dj » 1.25/? 


Experimental data (Powell 1953a) suggest that the shock-cell spacing is related to the feedback 
acoustic wave length as 

s!X~QA (11) 
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where X is an integer. 

Detailed measurements (Harper-Bourne 2002) suggest that the shock-cell spacing s is not 
uniform, but varies somewhat along the jet. A non-uniform shock-cell spacing, which may affect 
the directivity, is unaccounted for in the present theory. Experimental investigations of Massey et 
al. (1993) also indicate that the effect of heating on the directivity of screech is largely 
unimportant. No direct effect of jet temperature on the directivity of screech is however evident 
from Eq. (10). 

C. Existing Relations for Screech Frequency 

Tam et al. (1986) proposed for the screech frequency a semi-empirical expression of the form 

fdj/uj=f{Mj,r,To/T a ) (12a) 

where T a represents the ambient temperature. The fully expanded jet diameter dj is related to 
the jet exit diameter and the nozzle design Mach number M ^ as 

dj /d = f(Mj,M d ,y) (12b) 

This formula agrees well with the data and accommodates the effect of jet temperature (hot and 
cold jets). Massey and Ahuja (1997) provided separate expressions, similar in form to Eq. (12), 
for different screech modes depending on the Mach number range. 

The screech frequency increases with temperature since at a fixed Mach number, the 
convective velocity u c of the instability wave increases with temperature (Massey et al. 1993). 

The frequency is inversely proportional to the total time of travel in the screech cycle, (Tj + Ti ) , 

\ 

where 7j refers to the time for the downstream-propagating instability wave to traverse the 
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shock-cell system starting from the nozzle exit, and T 2 the time elapsed for the upstream- 
propagating feedback acoustic wave to reach the nozzle lip region. Thus, as the temperature 
increases, the value of T\ increases so that the frequency decreases. 

III. RESULTS AND COMPARISON 

In all the results presented here, except for the screech amplitude, all the sound pressure levels 
(dB) are taken realtive to 20 pi Pa. The quoted levels are in dB linear (unweighted). 

A. Streamwise Distribution of Acoustic Source Strength 

The streamwise distribution of the acoustic source pressure, as predicted by the present 
model for the normal shock case, is depicted in Fig. 6 for various Mach numbers. The successive 
amplification of the stream disturbances along the shock cells is evident, with the source strength 
increasing with the Mach number. A comparison of the predicted amplification rate with the data 
of Kaji and Nishijima (1996) for M j = 1 .64 and of Panda et al. (1997) for M j =1.19 and 1 .42 is 

also illustrated here. The data are deduced from the original measurements in order to provide a 
direct comparison. In this context, it would be instructive to briefly present the data in the 
original form. 

The streamwise amplitude measurements by Kaji and Nishijima (1996) for M y =1.64 suggest 

that the highest intensity is located in the shear layer of the third or the fourth shock cell, whereas 
inside the jet the acoustic intensity is low. Interestingly, the data reveal low-pressure dips 
(dimple spots) in the source strength. According to Kaji and Nishijima (1996) these dips are 
believed to be generated by the interaction between sources. There is a possibility that they could 
also be due to the expansion regions in the flow. 
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Detailed measurements by Panda et al. (1997) for a circular jet over a wide range of Mach 
number (Mj = 1.19 to 1.51) show similar trends for the axial distribution of source strength, 

with the peak values as a whole showing an increasing trend with M j , ultimately approaching a 

level of saturation followed by a relatively rapid decline in the screech levels (Fig. 7b). 

The foregoing data for the amplification rate are deduced here such that the relative source 
strength from the data at the first shock cell is taken as that given by the theory. Fig. 6 suggests 
that the proposed theory for the acoustic source strength M j = 1.2 is in close agreement with the 

data of Panda et al. (1997) for all the shock cells. However, at Mach numbers of 1.4 and 1.6, the 
theory overpredicts the amplification rates for downstream shock cells. This discrepancy is 
perhaps due to the assumption of normal shock. 

Numerical simulations by Manning and Lele (2000) for a single shear layer-compression 
wave interaction suggest that the radiated acoustic pressure amplitude approaches a saturation 
level at increased distances from the origin of the shear layer (Fig. 7c). The trends of the present 
model qualitatively agree with these results, although in the exponential region, the amplification 
rates appear to be comparatively high with respect to those in the present model. However, a 
direct quantitative comparison of this result with the present model appears difficult in view of 
the fact that a single compression wave is considered in the computation. 

B. Screech Amplitude 

Fig. 8 presents the variation of the excess screech amplitude as a function of Mj as provided 

by the foregoing theory. Experimental data reported by various investigators (including Kaji and 
Nishijima (1996), Massey and Ahuja (1997) Raman and Rice (1994), Raman (1997), Ponton & 
Seiner (1992), and Seto et al. (2006) are also displayed here in an effort to validate the physical 
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model. The theory for the normal shock case {(3\ =90 deg) suggests that the excess screech 
amplitude increases from about 17 dB at M y = 1.05 to about 46 dB at M } = 2.0, with the slope 
of the curve slowly diminishing with increasing Af y . The predictions for j3\ - 60 deg , addressing 

the sensitivity of the results to the shock wave angle, show about 7 dB reduction in the screech 
amplitude relative to the normal shock case. The theory seems to be in reasonable agreement 
with the data from various sources over a range of Mach number, with the normal shock 
providing the upper limit, and the oblique shock = 60 deg ) yielding the lower limit. 

Measurements for absolute screech amplitude (which includes broadband noise) for a 5:1 
aspect ratio rectangular nozzle obtained by Panda et al. (1997) for a wide span of M } confirm 

the general character of the present theory (Fig. 9). However, this report does not present the 
excess tonal noise for a direct (quantitative) comparison with the present calculations which 
provide only the excess screech magnitude. It is believed that the saturation (leveling) of the 
amplitude at high Mach number (say M > 4 ), as indicated by the data, may be qualitatively 
explained to some extent if the shock-wave angle varies (decreases) with downstream distance. 
The relatively rapid decline in the amplitude at M j = 1.7 is perhaps related to the change in the 

screech mode. 

The effect of screech mode on the screech intensity is beyond the scope of the present 
investigation. The dependence of screech intensity on jet temperature is also excluded from 
consideration in the present report. A detailed account of nozzle lip thickness is also not 
addressed here. 
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C. Screech Directivity 

Sample calculations of the directivity of screech for a typical value of Mj = 1 .4 for a three- 

source system are depicted in Figs. 10a and 10b, respectively, for the fundamental and the 
second harmonic. An effective value of M c = 0.7 is chosen as a representative value. Results for 

variable strength as obtained from Eq. (10) with c *2 =2.75, <33 =5.79 are compared with the 
case of uniform source strength (02 = <23 = 1 ). The values of <2 2 and <23 are obtained from Fig. 6 
at =1.4, relative to the first shock-cell strength of 3.5. For comparison purposes the 
directivity data of Norum (1983) for a choked nozzle at My = 1 .49 are included. The frequencies 

of the fundamental and the second harmonic in this data are respectively 19 Hz and 38 Hz. 

For the fundamental (Fig. 10a), the directivity pattern displays three lobes. One lobe, namely 
the principal lobe, radiates upstream and forms part of the feedback loop. The next dominant 
lobe radiates downstream at a relatively small angle from the downstream flow direction. The 
third lobe is the weakest one, which radiates in a direction normal to the jet axis. The directivity 
of the fundamental displaying the predominant lobe in the upstream direction is consistent with 
the well-known trends of the screech data (Norum 1983). It is evident that the effect of variable 
strength relative to uniform strength is displayed only in the weakest lobe and is not significant 
in the primary lobe directed upstream. The data for the fundamental however do not display the 
weakest lobe (normal to the jet axis). 

Fig. 10b for the directivity of the second harmonic suggests that the peak radiation due to 
screech is in a direction normal to the jet axis — an established result that is in concurrence with 
the available measurements, such as those of Norum (1983). Again, the variation of strength has 
some influence on the predicted directionality of the second harmonic in the three subordinate 
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lobes only and does not alter appreciably the directivity of the primary lobe directed normal to 
the jet axis. The theory considerably underpredicts the upstream lobe levels for the harmonic. 
Larger number of sources may in general improve the accuracy of the predictions for the 
directivity- both the fundamental and the second harmonic. 

The proposed method is applicable also to a forward flight stream surrounding the jet by 
appropriately modifying Eq. (2a) and (2b), and with a proper choice for the value of M c in Eq. 

(10). Measurements by Krothapalli et al. (1996) with a ffeestream Mach number from 0 to 0.32 
with an underexpanded heated supersonic jet (nozzle pressure ratios up to 4.5) reveal that the 
screech intensity is unaltered by the forward flight. Cold jet data by Norum and Shearin (1988) 
indicate the existence of intense screech tones in simulated forward flight up to Mach 0.41. 

IV. CONCLUSIONS 

Quantitative estimates for the excess screech tone amplitude over a range of fully expanded 
jet Mach number, as derived from the present model, seem to satisfactorily agree with the 
experimental data for underexpanded supersonic jets. The predicted streamwise amplification of 
the acoustic source strength is also found to be consistent with the available measurements, 
especially at low Mach numbers. This favorable comparison between the theory and the 
measurements suggests that the proposed physical model for the amplification of downstream- 
propagating disturbances of finite amplitude based on shock-acoustic wave interaction theory 
represents a plausible mechanism of screech in the feedback cycle. The results also indicate that 
the directivity of screech due to variable source strength of the discrete monopole source in each 
shock-cell does not appreciably differ from that based on uniform source strength insofar as the 
principal lobes are concerned. The proposed hypothesis for the determination of screech 
amplitude appears to hold promise in our effort to refine the jet screech modeling. 
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Captions to Figures 


FIG. 1. A typical narrowband farfield shock noise spectrum (adapted from Seiner 1984). 

FIG. 2. Schematic of the feedback loop in supersonic jet screech according to Powell (1953a). 

FIG. 3. Idealized distribution of acoustic sources in an underexpanded jet. 

FIG. 4. Interaction of an acoustic wave with a stationary shock wave (McKenzie & Westphal 
1958). 

FIG. 5. Acoustic pressure amplification as a function of upstream Mach number. 

FIG. 6. Streamwise amplification of acoustic source strength from present theory. 

FIG. 7a. Streamwise amplification of acoustic source strength from the data of Kaji & Nishijima 
(1996) for Mj =1.64. 

FIG. 7b. Streamwise amplification of acoustic source strength from the data of Panda (1997). 

FIG. 7c. Streamwise amplification of acoustic source strength from the numerical simulation of 
Manning & Lele (2000). 

FIG. 8. Comparison of the dependence of excess screech amplitude on Mach number. 

FIG. 9. Dependence of measured screech amplitude on Mach number, from data of Panda (1997). 
FIG. 10a. Directivity of the screech fundamental. 

FIG. 10b. Directivity of screech second harmonic. 
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FIG. 1 . A typical narrowband farfield shock noise spectrum (adapted from 
Seiner 1984). 
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FIG. 2. Schematic of the feedback loop in supersonic jet screech according to Powell (1953a- 
c) 



feedback 
acoustic wave 


shock-cell 


subsonic 


supersonic 

expansion compression 


stream dismrbance 


shear layer 



FIG. 3. Idealized distribution of acoustic sources in an underexpanded jet. 
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FIG. 4. Interaction of an acoustic wave with a stationary shock wave (McKenzie and 
Westphal 1958). 



FIG. 5. Acoustic pressure amplification as a function of upstream Mach number. 
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FIG. 6. Streamwise amplification of acoustic-source strength from present theory. 
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FIG. 7a. Streamwise amplification of acoustic source strength from the data of Kaji and 
Nishijima for A/ y =1.64. 
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FIG. 7b. Streamwise amplification of acoustic-source strength from the data of Panda (1997). 
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FIG. 7c. Streamwise amplification of acoustic-source strength from the numerical simulation 
of Manning and Lele (2000). 
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FIG. 8. Comparison of the dependence of excess screech amplitude on Mach 
number. 



FIG. 9. Dependence of measured screech amplitude on Mach number, from data of Panda 
(1997). 
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FIG. 10a. Directivity of the screech fundamental. 
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FIG. 10b. Directivity of the screech second harmonic. 
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